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DYNAMIC  PRODUCTION  NETWORKS 
by 

Ronald  W.  Shephard 

1.  INTRODUCTION 

Production  systems  are  typically  an  Interacting  collection  of  production 
activities  the  outputs  of  which  may  be  intermediate  products  serving  as 
inputs  for  some  activities  of  the  system  or  final  products,  or  both. 

The  flows  of  intermediate  products  may  be  represented  as  arcs  of  a  linear 
graph  with  activities  represented  by  nodes.  In  such  terms  production  may 
be  modeled  by  a  directed  network  connecting  activities  along  which  system 
exogenous  inputs,  intermediate  and  final  products  flow  dynamically.  The 
developments  to  follow  are  addressed  to  this  structure  of  a  dynamic  produc¬ 
tion  correspondence.  See  Shephard/FSre  [1980]  for  the  abstract  model  of  a 
dynamic  production  correspondence  relating  histories  of  exogenous  inputs 
to  histories  of  net  outputs  without  consideration  of  the  network  structure 
and  intermediate  products. 

In  the  development  of  the  abstract  structure  of  a  production  network 
as  a  dynamic  production  correspondence,  certain  network  axioms  are  needed 
to  characterize  the  role  of  intermediate  products,  and  to  verify  that  the 
correspondence  between  dynamic  flows  of  network  system  exogenous  inputs 
and  final  outputs  obeys  the  axioms  taken  for  such  systems  and  used  for  the 
network  activities  in  treating  intermediate  products  as  activity  exogenous 
inputs. 

With  this  theoretical  foundation  one  may  progressively  develop 
"ACTIVITY  ANALYSIS"  dynamic  models  from  abstract  to  computational  forms. 

The  computational  dynamics  for  production  networks  is  an  interesting  de¬ 
parture  from  the  superficial  practice  of  indexing  variables  by  time  and 


writing  down  the  Hamilton-Jacobi  differential  equations  for  minimizing 
or  maximizing  some  statement  for  an  unstructured  production  model.  It 
would  seem  that  representing  production  dynamics  economically  by  mechanical 
analogies  is  an  idle  practice,  despite  the  mathematical  artistry  of  the 
coordinate  and  variable  substitutions  involved  in  so  representing  the 
equations  of  motion  of  an  arbitrary  point  mass  system  subject  to  forces, 
inertia  and  equations  of  form  restricting  motion. 

In  the  discussion  of  computational  dynamics  for  production  networks, 
discrete  time  points  will  be  used  with  step  functions  for  histories. 

This  approach  enables  use  of  the  techniques  of  mathematical  programming 
for  practical  purposes. 

For  acyclic  production  networks  a  dynamic  computation  of  a  feasible 
output  trajectory  will  be  outlined  to  illustrate  the  evolutionary 
character  of  "motion"  for  production  networks,  and  to  show  possibilities 
for  Time  as  well  as  Factor  Substitution  in  production,  an  added  dimension 
for  dynamic  economics. 
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2.  ABSTRACT  STRUCTURE  OF  DYNAMIC  PRODUCTION  NETWORKS 


Consider  a  production  system  with  N  elementary  activities  and  a 
final  activity  (N  +  1)  for  recording  final  output  rate  histories.  All 
histories  are  taken  as  time  rates.  As  notation,  the  source  of  exogenous 
input  rate  histories  is  denoted  by  Aq  .  Thus  the  production  system 
consists  of  elementary  activities  A^.A^^ . ^’^+1  . 

The  primitive  elements  of  the  goods  and  services  related  to  this 
network  of  production  activities  are  time  rate  histories  of  system  exogenous 
inputs,  activity  outputs,  intermediate  product  transfers  and  final  outputs. 


each  defined  on  the  nonnegative  real  line  F+  .  Each  history  is  an 
element  of  (L<B)+  .  i.e.,  the  nonnegative  domain  of  equivalence  classes 
of  bounded  and  measurable  real  functions  defined  on  K+  .  (Two  functions 
are  equivalent  if  they  differ  only  on  a  subset  of  measure  zero.)  The  norm 
||f||  of  a  function  f  e  (Lw)+  is  taken  as  the  essential  supremum. 
Addition,  and  multiplication  of  histories  by  a  positive  real  number,  are 
taken  pointwise  in  time. 

As  notation: 


x  ■  (x-,x0,  ...»  x  )  e  (L  )°  is  a  vector  of  input  rate  histories 
for  n  system  exogenous  Inputs  (factors) . 

Vi  "  (Vi’Vi’  V’i)  e  (L»^+  »  "  1*2»  •••»  N)  are  vectors  of 

net  output  rate  histories  for  the  N  activities  of  the  network. 

ViJ  *  (Vij,Vij*  ****  ^j)  e  (L«»)+  »  1  “  L»2»  •••*  N  »  are  vectors 
of  transfer  rate  histories  from  A^  to  A^  as  intermediate  and 
final  outputs  for  j  *  1,2,  ...,  N  ,  and  j  ■  N  +  l  respectively. 

u  *  (u, ,u.,  . . . ,  u  )  e  (L  )®  is  a  vector  of  net  output  rate 

1  Z  01  00  T 

histories  for  the  production  system. 


A  common  space  LC^)®  is  taken  for  all  activity  output  rate  histories, 
transfers  of  the  same  and  network  output  rate  histories.  Not  all  component 
histories  of  ,  u  need  be  positive  on  subsets  of  R+  of 

positive  measure. 


11*11  *  M3*  I  l*J  ! 

i  1 
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By  defining  the  distance  p(f,g)  between  two  histories  to  be  the  essential 
supremum  of  |f  -  g|  ,  i.e.,  p(f,g)  -  ||(|f  -  g|)||  ,  the  model  of  produc¬ 
tion  is  defined  by  primitive  elements  in  metric  spaces.  Notions  of  closure, 
convergence  and  limits  for  these  elements  follow  naturally. 

The  activities  of  the  production  network  each  follow  net  output  and 
net  input  dynamic  production  correspondences  defined  respectively  by: 


«  x 


(L  ) 


m 


«'+ 


e  2 


2m 

+ 


(Vi,N+l’  Vij) 


e  (LJ 


2m 


Li(Vi,N+l’  Vij) 


(L 


where  XQi  ,  i  «  1,2,  . . . ,  N  ,  is  an  allocation  to  from  the  vector 

x  ,  subject  to 


i-1 


x  .  < 
oi  - 


and  denotes  a  summation  with  the  term  for  j  ■  i  omitted. 
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The  axioms  governing  the  activity  dynamic  correspondences  may  be 
found  in  Shephard/FMre  [1980]  and  will  not  be  repeated  here.  Since 
the  activities  of  the  production  network  are  primitive  production  elements 
for  the  system,  they  are  taken  to  relate  only  exogenous  input  rate  histories* 
to  net  output  rate  histories.  They  may  range  from  elementary  processes  to 
factories,  depending  upon  the  extent  of  aggregation  for  the  production 
system  studied. 

The  production  network  dynamic  output  correspondence  is  expressed  by: 


n  (LJT 

x  e  <LJ+  -  W*)  e  2 


*®(x)  -  jU  £  (L«>”  :  u  £  vi,^l‘  ^  Xoi  ^ 
(vi,K+r  vij) e  ffi(xoi'  vji)  • 


i  *  1, 2 ,  ...,N/. 


The  production  network  dynamic  input  correspondence  is  expressed  by: 

(L  )** 

u  c  <Lj“  -  ]LK(u)  -  {x  e  (Lj”  :  u  c  lPW(x)|  e  2  "  +  , 

M<“) '  jx  E  tt->2 ;  -oi-  vi,h+i  : u 

|xoi’  Vjlj  e  Li|Vt,N+l’  Vij  j  *  1  “  1,2 . N| 

In  these  expressions,  the  output  vector  of  each  activity  is  taken  net 
of  the  use  by  the  activity  of  its  own  products.  The  allocation  of  the 


'Here  the  input  rate  histories  exogenous  to  an  activity  may  span  both 
system  exogenous  inputs  and  intermediate  product  transfers. 


total  vector  x  to  activities  is  taken  freely  disposable,  without  loss 
of  generality.  Free  disposability  of  activity  final  products  is  not 
likely  to  be  generally  applicable,  but  is  taken  so  here  for  simplicity 
of  presentation.  See  Shephard  [1981  forthcoming]  for  such  complications. 
However,  input  and  output  histories  are  not  taken  freely  disposable  for 
the  activity  dynamic  production  correspondences.  More  on  this  in  the  next 


section. 


3.  NETWORK  AXIOMS 


It  is  important  to  investigate  whether  the  production  network  corre¬ 
spondences  1PK  and  ]LW  obey  the  properties  (axioms)  taken  for  the  dynamic 
production  correspondences  of  the  activities  of  the  network.  The 
properties  taken  for  3^  ,  do  not  propagate  for  3*1?  and  3.W  without 

some  properties  (axioms)  postulated  for  production  networks,  since  inter¬ 
mediate  products  are  involved.  The  network  axioms  are: 


1PW.1  For  each  component  history  u^  e  (La>)+  of  u  e  (L^^ 
which  is  not  entirely  intermediate  product,  there 
exists  a  vector  x  £  (L  )^  such  that  u  e  ]PlN(x)  with 

CO  -f. 

I|uil!  >  0  . 

ffW.2  For  each  activity  A.^  there  exists  a  subset  (i^i^  ....  ife}  , 

1  <  k  <  n  ,  of  the  n  exogenous  input  histories  xQi  which 

is  globally  essential  and  strong  limitational  for  the  output 

vector  V.  ,  i.e.,  x  .  has  to  be  essentially  positive, 
i  °LS 

j  *  1,2,  ....  k  ,  for  |  |vi|  |  >  0  ,  and,  if  |  |xoj.  .x^,  .... 

x  II  <  B  g  R  ,  there  exists  for  each  V°  a  positive  scalar 

oi,  “  "H-  1 

k 

0  depending  upon  B  and  V°  such  that  £  0V°  (0^  such 

that  (V±)j  j f  0^  (7°)^  .  j  -  1.2 . m)  for  all 

vi  £  l[  v)  • 


PK.2S  For  each  activity  A^  ,  {i^,^.  ...»  i^)  is  taken  globally 
totally  essential  and  strong  norm  limitational ,  i.e.,  xQi 

has  to  be  essentially  positive,  j  *  1,2,  ...»  k  ,  for  any 
subset  of  to  be  essentially  positive,  and  there  exists 

a  scalar  A  depending  on  B  and  such  that  | |v^| |  £  A 

when  | | xQl^ ,  ....  xol^l I  <  B  ,  for  all  . 
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PH.3  If  u  e  (L  )®  is  summable  in  each  component  on  E.  there 
exist  ,  i  =  1,2,  . ..,  N  ,  suramable  in  each  com- 


1PH.4  If  tu  :  ■  Max  (Ess  Sup  {t  e  1R+  :  u^(t)  >  0})  <  «  and 

x  ,  V  (i  -  1,2,  N  ,  j  -  1,2,  ....  N,N  + 1)  yield 

u  ,  then 


yi(t)  -  x±(t)  ,  t  e  [O.t^) 
y±(t)  -  0  ,  t  e  Uu,-h») 


(i  “  1,2,  « . . ,  n) 


w1j(t)  «  V±J(t)  »  c  e  t°»*u)  I1  *  1»2> 


wij (O  *  0  ,  t  e  Uu,+«) 


j  =  1,2,  ....  N.N  +  Ij 


PH. 5  For  each 


N 

activity  A±  ,  Vi>N+1  and  P  are 

weakly  disposable  subvectors  of  ,  i.e.,  if 

(w  j;  \) 8  j;  v) the" 


e1  e  [0,1]  ,  e2  e  [0,1]  . 


PH.6  For  each 


activity  At  ,  II  (xol,  £  Vyj  c  £  * 

(x\i-  “  j;  v3i) 8  Li(vi,mv  j;  vii)  *  > *  •  “ « ti 


N 


i.e.,  x  ^  and  }•'  V  are  weakly  disposable  vectors  of 
0  j«l  2 

input  rate  histories. 


A 


]PW.7  For  each 


activity  A±  if  ^VljN+1,  V±jj  e 


exists  for 


£  Vji)  3nd  !|Vi,N+l*  VijH  "  °  *  thete 
0  e  HR^  a  positive  scalar  Xg*^  such  that 

» •  (vm>  j;  vu) e  -ife1^ • 


1EW.1  For  all  vectors  V  -  (V, ,V_,  . ..,  V„)  such  that  u  is 

1  i.  N 

obtainable  from  the  related  production  network,  the 
efficient  subsets 

w  ■  j’  vji) :  (*oi-  j;  vji)  ‘  vv  • 

(y of  £  "ji)  1  W 

for  j;^))'2’ 


for 


are  uniformly  bounded. 


.m 


Concerning  the  first  of  these  network  axioms,  since  (L<JO)+  is  taken 
to  span  all  intermediate  and  final  output  histories,  some  component  histories 
of  u  may  be  essentially  null.  Hence  the  need  for  this  network  axiom. 

Axiom  1PW.2  states  that  the  activities  of  the  network  cannot  produce  net 
output  without  some  system  exogenous  inputs,  and  if  the  vector  of  essential 
system  exogenous  input  histories  is  bounded  in  the  norm,  the  scaling  of  each 
possible  vector  of  output  histories  is  bounded.  The  stronger  version  1PN.2S 


(2) 


The  sign  >_  means  that  at  least  one  component  is  "greater  than"  on  a 
subset  of  R+  of  positive  measure. 


replaces  global  essential  by  global  totally  essential,  and  strong  limita- 
tional  by  strong  norm  limitational.  See  Shephard/F£re  [1980]  for  definitions. 

Axiom  PM.  3  is  a  statement  for  the  network  correspondence  similar  to 
1L.T.1  taken  for  the  activity  correspondences. 

Axiom  PM. 4  is  a  statement  for  the  network  correspondence  like 
axiom  1L.T.2  assumed  for  the  activity  correspondences. 

Axiom  PM. 5  permits  scaling  of  activity  distributions  to  final  output 
independently  of  the  scaling  of  vectors  of  intermediate  product  transfer 
histories.  A  similar  assumption  is  made  in  PM. 6  for  vectors  of  exogenous 
input  histories  and  vectors  of  intermediate  product  input  histories  to 
an  activity. 

Axiom  PM. 7  postulates  that  the  scaling  of  a  vector  of  output  histories 
for  an  activity  may  be  obtained  by  the  same  scaling  of  intermediate  product 
input  histories  while  the  scaling  of  exogenous  inputs  may  exhibit  decreasing, 
constant  or  increasing  returns  to  scale. 

The  axioms  for  dynamic  production  correspondences  require  that  the 
subsets  be  bounded,  and  the  network  axiom  IN.l  results  in  this 

same  property  for  the  production  network  as  a  whole. 

Two  topologies  are  considered  for  the  metric  spaces  (L^)"  ,  (L^)”  , 

(L  )  ,  (L  )  .  The  norm  topology  under  the  essential  supremum  norm, 

and  a  weak*  topology  for  those  spaces  by  price  histories  taken  in  . 

With  the  foregoing  network  axioms,  the  axioms  for  the  activity 
dynamic  production  correspondences  propagate  for  the  network  correspondences 
PM  and  ]LN.  See  Shephard  [1981,  forthcoming].  The  closure  property  P.5 
is  the  only  property  with  complications  for  showing  propagation.  Two 


alternatives  arise: 


1P.2S,  IP. 5  (NORM  TOPOLOGY)  and  P.5  BIS  for  the  activities,  with 
PM.  2  and  PM.  5, 

P.2,  IP. 5  (WEAK*  TOPOLOGY)  and  P.5  BIS  for  the  activities,  with 
PM.  2  and  PM.  5. 

Regarding  the  production  network  efficient  subset 

PM(u)  *  |x  e  (L^)”  :  x  e  ILM(u)  ,  y  t  ILM(u)  for  y  <_  xj  , 

1N.1  implies  EM(u)  is  bounded.  1EM.1  is  required  because  there  are  an 
unbounded  number  of  ways  in  which  the  outputs  V  (i  *  1,2,  . . . ,  N) 
may  be  composed  to  yield  u  . 
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4.  AN  ABSTRACT  ACTIVITY  ANALYSIS  MODEL  FOR  DYNAMIC  PRODUCTION  NETWORKS 
In  some  physical  units  or  dimensionless  terms,  let  z  ■  ...» 

K 

zN)  e  (Lw)+  denote  a  vector  of  intensity  functions,  stating  for  each 

activity  the  intensity  of  operation.  These  intensity  functions  are 

-  N 

bounded  by  a  vector  z  e  (Laj)+  of  nonnegative  intensity-bound-functions 
expressing  the  Inherent  limitations  arising  from  physical  limitations 
not  otherwise  reflected  by  the  exogenous  service  Inputs  of  facilities, 
equipment,  and  also  by  product  design. 

Technical  coefficients  are: 


A  :  -  |  lA-'-A2  •••  AN|  I111  , 

*ij  I  ’  1  “  1,2» 


A1:  -  I  I a±1a12  ••• 


...  N  »  aij  G  (Loo)  + 


A  :  -  |  |ALA2  •••  AN||TO 
=i 


A  :  -  | ! ailai2  *•*  ®iml I  »  i  ■  1.2, 


..,  m  ,  a±j  e  (LJ+ 


t  *  -  || •••  , 


N.  i TR 


t  :m  I  I°ilci2  **’  Cimll  *  1  “  X’2’ 


■  •  ctj  e  (Lw)+ 


in  some  units  such  that  z^A*  ,  z^A*  and  z^0X  are  time  rate 
histories  of  exogenous  and  intermediate  product  inputs,  and  time  rate 
histories  of  outputs  respectively.  Then  the  production  network  dynamic 
output  and  input  correspondences  are: 


3PW(x)  ■  |u  e  (L^)”  :0<z<z,zA<x,  z(0  -  A)  >  0  , 

u  <  *($  -  A)| 

JLK(u)  ■  |x  c  (L^)”  :  0  <  z  <  z  ,  z($  -  A)  >  u  ,  x  >  za|  . 


LI 


4 


« 


ft 


» 


For  these  formulations  certain  assumptions  are  made  concerning  the 
technical  coefficients: 

For  each  activity  some  exogenous  input  is  required  during  [0,+“) 
except  on  subsets  of  measure  zero. 

Each  exogenous  input  is  required  by  some  activity  during  [0,+») 
on  subsets  of  positive  measure. 

Each  activity  can  produce  some  output  on  a  subset  of  [0,+")  of 
positive  measure. 

Each  output  is  produced  by  some  activity  on  a  subset  of  [0,-H») 
of  positive  measure. 

The  technical  coefficients  are  merely  taken  as  nonnegative 

as  stated. 

The  input  and  output  histories  have  been  taken  freely  disposable 
for  simplicity  of  expression.  Extensions  for  limited  disposability  will 
be  given  in  Shephard  [1981,  forthcoming]. 

As  abstract  statements  these  two  dynamic  production  correspondences 
do  not  provide  one  with  computational  systems,  but  point  the  way  toward 
such  dynamic  structures  for  production  theory,  expressing  the  inherent 
structure  of  production  dynamics. 


5.  A  COMPUTATIONAL  DYNAMICS  FOR  LEONTIEF-LIKE  PRODUCTION  NETWORKS 
Partition  the  nonnegative  real  line  into  half  open  intervals 

[t — l>t)  ,  t  ■  1)2,3,  • ■ •  . 

Take  the  intensity  functions  z  «  (z, ,z„ . z„)  constant  on  each 

1  L  N 

interval.  Thus  among  the  functions  (L  )  only  step  function  intensity 
histories  are  considered.  This  restriction  is  natural  for  a  computational 
system.  The  unit  of  time  is  arbitrary.  In  the  same  way,  histories  of 
system  exogenous  inputs  available  will  be  taken  as  step  functions. 

At  this  juncture  it  is  convenient  to  introduce  more  detail  than 
previously  considered.  Exogenous  inputs  may  be  storable  as  well  as  non- 
storable,  but  not  both.  In  order  to  accommodate  this  fact  the  notation 
for  exogenous  input  histories  is  expanded  to 

x  C  (L  )S  STORAGE 

oo 

y  e  (L 00)n”S  NONSTORAGE. 

Initial  inventories  of  storable  exogenous  inputs  are  included  in  the 
system  exogenous  input  histories 

x^(0)  ,  i  "  1,2,  ...»  s  • 

Concerning  Intermediate  products.  Initial  Inventories  of  the  same 
are  denoted  by 

v®  ,  k  ■  1,2,  ...,  m 


shared  by  all  activities  of  the  network,  and  capacities  for  storage  of 
intermediate  products  are  denoted  by 
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ok(t)  k  m  1,2,  ■  •  • ,  in  ,  t  *  1,2,3,  ••• 


independent  of  activity  producing.  Intermediate  and  final  outputs  for 
production  during  [t-l,t)  are  taken  to  be  realized  at  t  .  The 
intensity  functions  z^(t)  and  intermediate  product  inputs  are  applied 
at  (t  -  1)  for  [(t  -  l),t)  . 

Since  some  activities  may  yield  outputs  which  are  both  intermediate 
product  and  final  product,  the  intensity  functions  z  e  (L^)^  are  bi¬ 
furcated  as 


Iz^  e  (Los)+  »  intermediate  product  production 
Fz^  e  (Loo)+  ,  final  product  production. 


In  these  terms  the  computational  dynamics  of  the  production  network 
is  subject  to  the  following  constraints:  (See  Shephard  et  al.  [1977]) 

(1)  zt(t)  -  IZl(t)  +  Fz±(t)  ,  zt(t)  <  ^(t)  , 

Iz±(t)  >  0  ,  Fzt(t)  >  0  ,  (i  -  1,2,  ...,  N)  ,  (t  -  0,1,2,  ...) 


(2) 


t  N 


l  l  a  (t)z  (t)  <  l  x  (t)  ,  (j  -  1,2,  ...,  s)  ,  (t  -  0,1,2,  . 
T«0  1*1  3  1  T-0  3 

N 

(3)  l  a. .  (t)z.  (t)  <  y  (t)  ,  (J  -  s+l,s  +  2 . n)  ,  (t  -  0,1,2,  .. 

i-1  3  3 


N  t  N  t-1 

(4)  l  t  m  vk  +  E  I  c  (x  +  Dlz  (T) 


Z  aik(0)zi(0)  <  v°  ,  (k  -  1,2 . m)  ,  (t  -  1,2,3,  ...) 

N  t-1  N  t 

(5)  v?  +  l  l  c  (t  +  1)Iz  <t)  -  l  l  »ik(T)zi<T>  i  °k(t>  » 

*  i-1  T-0  1  i-1  T-0 

N 

Vk  “  ^  *ik(0)*i(0)  -  °k<0)  ’  <k  “  1»2»  •••»“>  »  *  1»2»3» 


•  •  • 


The  constraints  (1)  limit  the  intensities  of  operation  of  the 
activities  of  the  production  network  to  be  nonnegative  and  not  to  exceed 
certain  physical  limitations  not  otherwise  expressed  by  the  service 
components  of  the  exogenous  input  histories. 

The  constraints  (2)  and  (3)  limit  exogenous  inputs  to  available 

3  n_8 

resources  as  given  by  the  time  histories  x  e  (L01l)+  and  y  e  (1.^)+ 

The  constraints  (4)  require  that  inputs  of  intermediate  products  by  the 
network  activities  do  not  exceed  the  supply  available  from  the  outputs 
of  the  activities  and  initial  inventories.  The  last  set  of  constraints 
does  not  allow  accumulation  of  inventories  of  intermediate  products 
beyond  capacities  for  the  same.  All  these  constraints  are  limitations  on 
network  activity  intensities. 

The  foregoing  constraint  system  does  not  in  any  way  predetermine 
the  allocation  of  resources.  However,  it  is  a  general  basis  for  a  com¬ 
putational  dynamics  of  Leontief-like  production  networks.  But  since  nothing 
has  been  specified  concerning  final  outputs  the  dynamic  evolution  of  the 
production  network  is  not  directed.  The  intensity  functions  Fz^  , 

(i  *  1,2,  ...,  N)  are  so  to  speak  free  elements. 

Since  the  output  set  ]PK(x,y)  in  general  may  be  taken  to  exhibit 
weak  disposability  for  vectors  u  e  (L^)™  of  output  histories,  one  may 
generate  for  various  feasible  output  mixes  of  output  histories  the  maximal 
scalar  extension  in  lPK(x,y)  ,  i.e.,  determine  the  dynamic  evolution  of 
the  system  to  find  points  on  the  frontier  of  2Pli(x,y)  .  Also  for  some 
price  histories  of  the  various  outputs  and  inputs,  one  may  seek  to  control 
the  dynamics  of  the  system  to  maximize  (Revenue-Cost  of  Resources) .  In 
either  case  it  is  convenient  to  make  these  generations  of  the  dynamics 
of  the  system  with  respect  to  a  finite  planning  horizon  T  . 


Now  suppose  there  are  1  <  P  <  m  net  products  possible  for  the 


system.  Let 

*N+l,p  £  K+ 

be  coefficients  related  to  an  index 


z  e 
N+l 


such  that 


*N+l,r  *3*1  ’  P  “  1,2 . p 

defines  the  total  amount  of  the  pth  final  product  accumulated  over  the 
planning  interval  [0,T]  .  The  coefficients  a^  p  may  be  chosen  to 
determine  a  specified  output  mix.  Then,  in  order  to  direct  the  system 
this  way  one  adds  the  constraint 

(6)  Vl.p  *  ZN+1  ^  cip(t)FZi(t-l)  ,  (P  -  1,2,  ....  P) 

and  determines  the  dynamic  evolution  of  the  system  by  the  following 
linear  program: 

(  MaX  ZN+1 

| Subject  to  2^+1  *  0  »  (1)»(2),(3),  ...»  (5), (6)  on  [0,T]  . 

A  calculation  of  this  type  has  been  outlined  by  Leachman  [1980]  for 
the  case  of  a  network  with  each  activity  producing  a  single  output,  possibly 
with  cycles  in  the  network,  i.e.,  for  a  Leontief-like  network.  Also  similar 
calculations  were  outlined  for  this  particular  network  to  minimize  the  cost 
of  obtaining  a  given  program  of  output  histories. 


It  Is  clear  from  the  foregoing  that  productivity  of  a  production  system 
is  a  complicated  concept.  For  any  given  output  mix  one  may  determine  the 
maximal  throughput  for  the  mix,  but  this  will  vary  with  output  mix.  The 
resource  inputs  likewise  condition  the  throughput.  How  is  one  then  to 
determine  properly  the  maximal  productivity  potential  of  a  system  unless 
resources  and  output  mix  are  balanced  for  this  purpose?  Obviously  a  good 
deal  of  study  is  needed  for  such  questions.  With  declining  energy  resources 
it  would  be  of  interest  to  know  how  the  maximal  throughput  of  various  output 
mixes  would  retrogress. 

In  the  general  terms  of  the  formulation  of  production  networks  used 
here,  real  capital  is  expressed  in  terms  of  service  input  rates,  and  for 
those  components  of  (x,y)  which  are  fully  utilized  in  the  dynamic  solution 
of  the  linear  program  and  reflect  real  capital,  one  may  make  marginal  analyzes 
(linear  programs)  for  incremental  changes  in  such  inputs.  Indeed  for  energy 
supply  decreases  one  may  seek  to  estimate  the  consequences  for  economic 
sectors  by  the  dynamic  Leontief  model  detailed  by  Leachman,  using  "constant 
dollar  values"  for  measures  of  aggregate  output. 

In  order  to  drive  the  dynamics  of  the  system  by  maximizing  total  revenue 
minus  cost  of  resources  (x,y)  applied,  one  need  merely  define  for  given 
price  histories  r ^(t)  »  qj(t)  the  objective  function 


n(T) 


N  /  T-l  \ 

l  I  l  r  (t+l)c  (t  +  l)Fz(t)| 
i-l\t-0  p  ip  / 

-1(1  q.<T)x  (t)  -  l  q,(r)y  (x)J 
t-0  \j«l  3  3  j-s+1  3  3  / 


and  determine  the  dynamic  evolution  of  the  system  by  the  linear  program: 
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A _ t 


{Max  n(T) 

Subject  to  (1)  ,  (2)  ,  (3) . (5)  on  [0,T]  . 

Here  x^ (t)  ,  y (t)  are  also  variables  in  the  linear  program.  Obviously 
refinements  can  be  made  in  this  kind  of  generation  of  the  dynamics  of  the 
production  network  by  discounting  and  adding  other  costs  to  complicate 
the  definition  of  profit.  For  the  purposes  of  this  paper,  interest  is 
mainly  in  structure,  and  not  management  practice. 


L.  k 


j 
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6.  ACYCLIC  PRODUCTION  NETWORKS 

In  a  production  network  without  cycles  the  activities  of  the  network 
may  be  ordered  so  that  the  outputs  of  an  activity  can  serve  as  inter¬ 

mediate  product  only  for  Ai+i»Ai+2*  ***,  \  »  or  w*iat  the  same, 
intermediate  product  Inputs  to  A^  can  arise  only  from  A^.Aj,  A^_^  . 

Then,  for  preassignment  of  (1)  resources,  (2)  preallocation  of  distributing 
intermediate  products  and  (3)  preallocation  of  activity  outputs  to  final 
product,  dynamic  trajectories  may  be  calculated  (evolved)  for  the  intensity 
functions  of  the  activities  and  thereby  determine  the  time  histories  of 
intermediate  and  final  products.  See  Shephard/et  al.  [1977]  for  such  cal¬ 
culations  for  shipbuilding. 

For  the  general  model  of  computational  dynamics  for  production  networks, 

let 


*01  ■  <*oir*oi2 . W  1  *  1*2 . " 

yol  ’  <yoll,yol2 . yoi(n-.))  1  "  1,2 . * 

denote  a  preallocation  of  exogenous  input  histories  to  the  activities,  with 


N  N 

7  x  <  x  ,  T  y  .  <  y  . 
i-i  01  ■  iii  oi  ■ 


Further,  the  ordering  taken  for  the  network  is  such  that  takes  the 

form 


Isu5i2  si,i-i  •  0  “•  °M  1  “  ly2 . N  ■ 


Let 

,  (j  -  1,2 . H)  ,  (k  -  1,2,3 . m)  ,  (t  -  1,2,3,  ...) 

J  f  * 
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t 


denote  the  fraction  of  the  output  rate  history  of  going  to  net 

output  at  the  time  t  .  These  preassignments  satisfy 


M4.I 

0  -  6j  k(t)  «  1  *  (t  "  1*2*3’  •••>  »  <k  “  1«2 . m)  ,  (j  =  1,2 . N) 


N 


j-1 


N+l 

5  j*(t)  >  0  ,  for  some  k  e  {1,2,  ...,  m}  ,  t  e  {1,2,3,  ...} 

Ji* 


and  the  net  output  histories  u  e  (L^)™  are  given  by 


N 


u,(t)  "  I  ifJWZ.U-Dc,  .  (t)  ,  (t  -  1,2,3,  ...)  ,  (k  -  1,2,  ....  m) 
j*l  2  ■3’ 


Let 


1  -  Aijk(t)  «  °  *  <k*1>2*  •••»«)  >  (j  “1*2,  . . . ,  N  -  1)  ,  (i-  (j  +l),(j  +2) 

denote  a  preallocation  of  the  k  output  rate  history  of  at  time  t 

the  activities  A^  ,  (i  =  (j+l),(j+2),  ...,  N)  as  intermediate  product. 
These  coefficients  satisfy 


I  A  (t)  -  1  -  5?T(t)  ,  (k  -  1,2,  ....  m)  ,  (t  -  1,2,3,  ...) 
i-j+1  1JK 

(j  -  1,2 . (N  -  1)  )  . 


Since  intermediate  product  time  histories  are  preallocated,  initial 
inventories  and  capacities  for  storing  intermediate  products  are  now  dis¬ 
tinguished  by  activity  as: 


v”k  (i  -  1,2,  ....  N)  ,  (k  -  1,2 . m) 


oik(t)  (i  ”  1*2,  ...»  N)  ,  (k  ■  1,2,  ...,  m)  ,  (t  ■  0,1, 2, 3,  •».)  . 


,  • .  • ,  N) 

to 


Note  that 


(0)  . 


Certain  simplifying  assumptions  may  be  taken  for  the  coefficients 
*ik(x)  *  cjk(T)  : 


either 

a.,  (t)  >  0  for  all 

t  e  [0,®) 

or  a  (t)  =  0  for 

all  t 

e  [0,")  ,  (i  -  1,2,  . 

...  20 

either 

c  (t)  >  0  for  all 

T  £  [0,+") 

or  c  ,  (t)  =  0  for 
J  K 

all  t 

e  [0,4«)  ,  (j  -  1,2, 

....  N)  . 

In  other  words,  an  activity  is  consistent  in  the  use  of  intermediate 
products  and  production  of  outputs.  This  does  not  preclude  alternative 
activities  with  different  intermediate  product  inputs  and  different  output 
commodities.  If  an  activity  cannot  satisfy  the  two  assumptions  it  may  be 
subdivided  until  it  does. 

Also,  with  little  if  any  loss  of  generality  one  may  assume 
o^j^(t)  >  n^^(t — 1)  ,  (t  *  1,2,  •  ••)  ,  (i  ®  1,2,  ...,  N)  ,  (k  *  1,2,  ...,  m) 


By  definition 

Aijk(T)  *  0  if  iik(T)  "  0  * 

One  does  not  transfer  intermediate  product  output  to  an  activity  not  using 
it. 

Then  the  previous  constraints  for  the  dynamic  system  take  the  following 


form: 
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(1)  0  <  2i(t)ii(t)  ,  (t  -  0,1. 2, 3,  ...)  ,  (i  -  1,2,  ....  N) 

*  t  t 

(2)  I  ,aik^T)zi^x)  -  £  Xoik^T)  ’  (k  =  1,2 . s>  ’ 

T“0  T»0 

(i  -  1,2,  ....  N)  ,  (t  -  0,1,2,  ...) 

(3) *  aik(t)z1(t)  <  yQik(t)  ,  (k  -  (s  +  1) , (s  +  2) ,  ....  n)  , 

(i  -  1,2,  ....  N)  ,  (t  «  0,1,2,  ...) 

<4>*  Jo  iik(T)zi(T)  =  Vik  +  Tl±  ^  AiJk(T)Cj,k(T)zJ(T’1) 

aik(0)zi(0)  -  vik  ’  "  1’2»3’  •••>  »  (i  "  2»3’4 . N)  » 

(k  -  1,2,3,  ....  m) 

(5)*  Vlk  +  j/  Aljk(T)cjk(T)zJ(T"1)  ‘  J0  aik(T)zi(T)  ^  °ik(T) 
(i  -  1,2,  ....  N)  ,  <k  -  1,2 . m)  ,  (t  -  1,2,3,  ...)  . 


Notice  now  that  the  intensity  functions  need  not  be  separated  into 

N+l 

Iz  ,  Fz  due  to  the  coefficients  6  ,  (t)  ,  (t  =  1,2,  ...)  . 

J  »*• 

It  is  of  some  interest  to  consider  the  dynamic  development  of  the 

trajectories  for  final  outputs  from  given  preallocations.  The  i»w<t'«i4iity 
*  *  * 

system  (1)  , (2)  ,  . . . ,  (5)  may  be  used  to  develop  trajectories  for  the 
intensity  functions  z^(t)  ,  t  -  0,1,2,  ...  which  generates  the  histories 
of  all  outputs,  both  intermediate  and  final.  A  policy  of  applying  feasible 
intensities  z^(t)  is  needed  in  order  to  get  a  specific  trajectory.  It 
is  convenient  to  seek  to  take  this  policy  as:  "utilize  the  maximal  value 
of  z^t)  possible  for  each  t  ■  0,1,2,  ...,"  i.e.,  a  Greedy  Policy. 

This  policy  ignores  the  problem  of  "variable  loading."  However,  it  will 
enable  one  to  observe  that  dynamic  systems  of  production  involve  Time 
Substitution  as  an  additional  dimension  beyond  those  for  static  or  steady 


state  systems. 
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The  Greedy  policy  for  this  system  is  generated  as  follows. 
Define: 


S(±)(t)  •  {k  e  {1,2 . s}  :  aik(t)  >  0} 

S^(t)  ■  (k  e  {s  +  l,s  +  2,  ...»  n}  :  >  0} 

E^(t)  -  {k  e  {1,2,  ...,  m}  :  aifc(t)  >  0}  • 


Define: 


R<« 


(0) 


if  S(i)(0)  +  0 
if  S(i) (0)  *  0 


R(±) (t) 


keS 

m  +00 


Min 

(i) 


T-0 


oik 


t-1 

(t)  -  l  a  (t)z  (t) 
-0  1 


(t) 


aik(t> 


if  S(i)(t)  *  0 
if  S(i)(t)  -  0 


R(1> 


(t) 


if  S(i)(t)  i  0 
if  S(i)(t)  -  0 


if  I(1)(t)  + 
if  I(i)(t)  - 
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Define 


QA(t)  “  Min  {R(i)(t),R(i)(t),W(1)(t),Zi(t)}  . 


Then  the  Greedy  Policy  is 


ZA(t)  -  Q1(t)  (i  -  1,2,  ....  N)  ,  (t  *  1,2,3,  ...) 

*  * 

if  constraint  (5)  is  satisfied.  However,  (5)  may  not  be  satisfied  at 
some  time  t  >  0  for  one  or  more  commodities  k  and  activity  i  .  Let 
k  yield  the  maximal  value 


(i-1) 

l  (t)c..  (t)*.(t-1) 

j=l  3o  1 o  3 


l  aik  (T)2i(T)  ‘  °ik  (to> 

T»0  O  O 


>  0 


Presumably  (5)  is  satisfied  at  (tQ  - 1)  for  i  and  kQ  ,  which  implies 


(i-1) 


l  A  .  (t  )c  ,  (t  )z  (t  -  1)  -  a.,  (t  )z  (t  )  >  0  . 
,L.  ilk  o  Ik  o  1  o  ik  o  i  o 

j“l  J  o  J  J  o 


Evidently,  (5)  may  be  satisfied  at  tQ  if  z^ (t)  *  0  for  t  *  (tQ  - 1)  , 

(j  -  1,2,  ...,  (i-1))  .  As  an  approximation  to  a  Greedy  Policy  one  may 
take  these  values  for  z^(tQ-l)  for  (j  ■  1,2,  ...,  (i-1))  . 

The  solutions  for  z^(t)  ,  (j  “1,2,  ...,  (i-1)  ,  (x  *  tQ  ,  (tQ  + 1) ,  (tQ  +  2) , 
...)  need  to  be  updated.  Then  one  may  proceed  t  with  the  calculation  for 
z±(to)  ,  zjL(to  +  l)  ,  etc. 

The  inventory  capacity  bounds  o^(t)  may  require  considerable  re¬ 
calculation.  However,  the  routine  suggested  is  a  simple  one  and  does 
provide  a  ready  policy. 


The  evolutionary  character  of  the  trajectory  for  a  greedy  type 


solution  is  evident.  In  the  order  (i  *  1,2,3,  ...,  N)  one  calculates 
z^(t)  for  (t  ”  0,1, 2, 3,  ...)  corresponding  to  the  given  preassignments 
of  activity  outputs  to  other  activities  and  final  output.  The  net  output 
histories  for  the  production  network  are  given  by 


^(t) 


N 

I 

j-1 


4jT1><t)vt),j(t-1> 


c 


1,2, 

1,2,3 


Now,  insofar  as  the  modified  Greedy  Policy  implies  for  some  activity 
intensities  that  periods  of  zero  intensity  occur,  one  may  seek  to  smooth 
load  such  production  intensities  by  operating  them  at  less  than  maximal 
intensity  without  altering  the  output  histories  obtained  over  some  planning 
period,  i.e.,  there  may  be  possibilities  for  Time  Substitution. 

Time  substitution  for  the  activity  intensity  trajectories,  and  through 
them  the  application  of  exogenous  input  and  intermediate  product  histories 
is  an  added  dimension  for  dynamic  production  networks  over  that  of  static 
or  steady  state  models  of  production.  The  related  problems  of  smooth 
loading  are  the  gist  of  dynamic  production  planning.  See  Leachman  [1980] 
for  some  mathematical  programming  methods  to  smooth  load  for  shipbuilding. 

Maximizing  throughput  as  discussed  in  the  previous  section  and  the 
one  to  follow,  and  other  optimizations,  are  not  always  meaningful,  because 
preassignments  of  exogenous  inputs  and  intermediate  and  final  outputs  may 
be  required  for  the  facilities  (activities)  of  a  large  production  network. 
Also  a  compounding  of  the  complexity  of  an  optimization  may  be  involved 
to  assure  a  reasonable  smooth  loading.  There  is  some  advantage  to  simulate 
a  dynamic  trajectory,  in  the  face  of  uncertainty  of  Information,  for  policy 
formation  which  can  be  updated  from  time  to  time. 


* 


A  COMPUTATIONAL  DYNAMICS  FOR  GENERALIZED  PRODUCTION  NETWORKS 


In  Sections  4,  5  and  6  above,  the  activities  of  the  network  were 
driven  by  intensity  functions  z^  e  (Los)+  operating  on  fixed  technical 
coefficients  for  system  exogenous  inputs,  intermediate  product  inputs  and 
activity  outputs,  independently  of  the  intended  distribution  of  activity 
outputs . 

Now  instead  of  a  single  intensity  function  z^  e  (I'00)+  for  A^  , 

let 

,  .  ,T  .N+1 

zi  :  "  zil,zi2*  ziN,Zi,N+l)  e  (IW+ 

denote  an  intensity  vector  with  z e  (La>)+  denoting  the  intensity  of 
operating  A^  for  output  to  go  to  A  ,  (j  ■  1,2,  ....  N,N  +  1)  .  The 
system  exogenous  input  histories  implied  by  z±  are  where 

A±:  =  TRANSPOSE  |  | •••  I  I  .  C1  =  •••»  N> 

Aij  5  *  ^aijlaij2  *“  aijn^  »  (j  "  1»2»  •••»  N’N  +  1)  • 
Similarly,  intermediate  product  inputs  are  given  by  z^A^  where 

A£  :  ■  TRANSPOSE  |  IS^A^  •••  AjA  (**].)  I  I  »  *  1>2»  N> 

^ij  Haijlaij2  ***  aijn^  »  “  1,2 . N.N  +  l) 

and  activity  outputs  are  given  by  z^0^  (i  ■  1,2,  ...,  N)  where 

t±  :  -  TRANSPOSE  || ***  Wi(N+l)'  ^  *  (i  "  1»2»  •••»  N> 
im  !fCijl’Cij2  •'*  ^jn11  »  (J  "  X’2 . N’N  +  1>  * 
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The  constraints  for  the  production  network  are: 


(1) 


** 


<♦>**  l  l  <  v°  +  l  l  c  (t)z  (t-1) 
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The  "motion"  of  this  system  for  a  production  plan  needs  direction  by 
optimizing  some  objective.  Maximum  throughput  of  a  given  output  mix  over 
a  planning  period  T  can  be  used  to  direct  this  production  system  by 
the  following  linear  program 


MaX  ZN+1 

Subject  to:  ^  >  0  ,  (1)**,  (2)**,  (3)**,  (4)**,  (5)**,  (6)**  on  [0,T] 


where 


(6) 


**  . 


N  T 


*N+l,kZN+l  -  Ci,N+l,k(T)zi,N+l(T  _1>  *  (k  “  1,2»  •**’  P) 


The  conventions  of  Section  5  are  used  in  this  formulation. 


The  minimal  time  horizon  for  a  production  plan  to  yield  given 
accumulations 

i  (k  *  1,2,  • • • i  P) 

** 

of  the  final  products  can  be  found  by  replacing  (6)  by 


(7) 


N  T 

Jl  Ci,N+l,k^T>Zi,N+l<T  ”  1} 


(k  -  1,2, 


P) 


**  **  **  ** 

and  consecutively  seeking  a  feasible  solution  to  (1)  ,  (2)  ,  (3)  ,  (4)  , 

** 

(5)  ,  (7)  for  increasing  T  *  1,2,3,  ...  until  feasibility  is  first 

attained.  See  Leachman  [1980]  for  this  computation  for  production  networks 
with  activities  each  of  which  produce  a  single  output. 

Other  optimizations  are  possible  to  drive  the  "motion"  of  the  production 
system;  and  additional  constraints  may  be  significant.  For  example,  large 
changes  of  service  inputs  may  involve  learning  and  shift  the  system  retro- 
gressively  to  higher  amounts  required  per  unit  intensity. 
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